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Prediction of Heated Jets 


e Jets with significant temperature differences have 
many important applications 
¢ Aeroacoustics 
° Cooling flows 
¢ Fuel injectors 
¢ IR signatures 


¢ Standard CFD methods (RANS) do a very poor job 
predicting these flows 
¢ Possible reasons 
¢ Turbulence model 


¢ Turbulent Prandtl number variation 
¢ Turbulent heat flux modeling 
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Turbulence Closure 


¢ Reynolds stress 
¢ Boussinesq approximation 


¢ Turbulent heat flux 
¢ Gradient diffusion & Reynolds Analogy 
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Approach 


¢ Use large-eddy simulation (LES) to examine the 
turbulent heat flux vector, q',, and turbulent 
Prandtl number,Pr,, in hot jets 


¢ Validate the LES using experimental data 
¢ PIV data for velocity 
¢ Rayleigh scattering and Raman spectroscopy for 
temperature 
¢ Make a leap of faith that if u’ and T’ are validated 
separately, then u’T’ should not be too bad 


¢ Compare results to Reynolds-Averaged Navier- 
Stokes (RANS) simulations and evaluate 


Large-Eddy Simulations 


¢ WRLES code 
¢ Explicit high-resolution finite-difference code 


e 11-pt DRP differencing scheme (Bogey & Bailly, 2004) 
with matching filter 


¢ 4-stage, 3 order Runge-Kutta time stepping 
¢ Hybrid MPI/OpenMP parallelization 
° Grid 
e Structured grid 
¢ 36 million points 
¢ 912x184x181 points downstream of nozzle exit 


RANS Simulations 


¢ Wind-US 
¢ Finite-volume 
e Structured grid axisymmetric mode 
¢ 2"¢_order upwind biased RHS 
¢ Full block-implicit LHS 
¢ SST-V turbulence model (vorticity based production 
term) 
° Grid 
¢ Taken from turbmodels.larc.nasa.gov 
¢ 73,151 points 
¢ Downstream of the nozzle exit 257x251 points 
¢ Provides grid converged solutions with Wind-US 


Round Jet Experiments 


¢ Small Hot Jet Acoustic Rig (SHJAR) 
e 2-inch nozzles: ARN2 and SMCOO0O 


¢ PIV Velocity Data 
¢ Bridges and Wernet, NASA TM 2011- 216807 
¢ Concensus dataset 
¢ Verified against hotwire and LDV 


¢ Rayleigh Scattering Temperature Data 
e Mielke et al, AIAA Journal, Vol. 47, No. 4, 2009 
¢ Point measurement SMCO00 


¢ Raman Spectroscopy Temperature Data 
e Locke and Wernet, NASA TM 2017-219504 
¢ Point measurement 


LES Methodology 


¢ Implicit LES 

¢ Nozzle boundary layer 
¢ No attempt to resolve a turbulent boundary layer 
¢ Transition occurs quickly in mixing layer 

¢ Non-dimensional time 
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¢ Startup time: 60t° 

¢ Averaging time: > 180t° 
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Flow Conditions 
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Axial Turbulence Intensity 
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Radial Turbulence Intensity 
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Radial Profiles —u’ & v’ 
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RMS Temperature 
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Radial Profiles — T’ 


3.0 

2.5 
SP 3 

2.0 
M, = 0.5 as 
T/T.. = 0.950 “1, 

0.5 

0.0 

3.0 

oT x/Dj = 2.0 x/Dj = 4.0 x/Dj = 6.0 x/Dj = 8.0 x/Dj = 10.0 x/Dj = 12.0 
SP 23 20 oO ao 

O° Rayleigh 

M, =0.5 


ti.=t77 


SP 27 af Pees 
M, = 0.9 O° Rayleigh 
Via i7e4 


SP'3 


M,=0.5 
T,/T.. = 0.950 


SP 23 


M,=0.5 
T,/T.. = 1.764 


SP 27 


M, = 0.9 
1/To= 1764 


x/Dj = 10.0 


x/Dj = 12.0 


0.05 
qi(Rj Uj dTj) 


x/Dj = 10.0 x/Dj = 12.0 


<<“ 
qi/(Rj Uj dTj) 


x/Dj = 10.0 x/Dj = 12.0 


Contours of <guu> and q, pes 


val 


+ Wik 
Se a 2) 
0.030 
0.024 
0.018 
0.012 
0.006 
0.000 


0 _8 
x/Dj 
~ AT Sane i 
q,, A UA | i] =<Ou | - 
3 i 
1 ya 
25 0.020 
0.015 
0.010 
1 &j 0.005 
0.000 


Contours of <puv> and q, pes 


V 


< 


<ruv>/ Ri UjA2 
a une) 


0.009 
0.006 
0.003 
0.000 


3&5 
/ Ri Uj dTi 
(A OST) 
0.012 
0.009 
0.006 
0.003 
0.000 
0 0 . ‘ 


_8 
x/Dj 


Turbulent Heat Flux 


¢ Turbulent heat flux model 
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¢ Radial component 
¢ RANS & LES agree surprising well 
¢ Mean temperature gradient is in radial direction 


¢ Axial component 
e LES predicts heat flux larger than radial component 
¢ RANS model predicts almost no heat flux 
¢ No temperature gradient in this direction 


¢ LES heat flux agrees with experiments in the literature 
¢ Magnitude 
* Fabris (1979): <uT> & <vT> similar in magnitude 
¢ Tavoularis & Corrsin (1981): <uT> larger than <vT> 
¢ Alignment (angle between temp. gradient and heat flux vector) 
* Current LES: 57° 
* Tavoularis & Corrsin (1981): 63° 


¢ Gradient diffusion model is not appropriate for this flow 
¢ Heat flux behavior is analogous to momentum flux 


Effect on the Energy Equation 


¢ Energy Equation 
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¢ Quantify the contribution of the missing axial 
component 
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Turbulent Prandtl Number 


¢ Treated as a constant but varies, 0.5 < Pr,< 1.0 
¢ Pr, = 0.7 is standard value for jets 


¢ Variable Pr, models often cited as a a solution to 
these types of problems 

¢ Yoder’s (2016) recent results showed no advantage 
for jets 


¢ Can be computed from the LES 
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Summary and Conclusion 


¢ LES and RANS methods were used to compute 
heated jet flows 


¢ RANS under-predicts spreading rate and inviscid core 
length (expected result) 


¢ LES agrees well with experimental data 


¢ Turbulent heat flux 
e LES results consistent with literature 
¢ RANS model fails to replicate physics 
¢ Gradient diffusion assumption not appropriate for jets 


¢ Turbulent Prandtl number 
¢ Little variation within the jet mixing layer 
¢ Pr, = 0.7 is consistent with literature 


